Relativistic hydrodynamics from 
Boltzmann equation with modified 
collision term 

Amaresh Jaiswafl, Rajeev S. Bhalerao and Subrata Pal 

£C) ' Tata Institute of Fundamental Research, Mumbai, India 



o 

(N 



Abstract 

Generalizing the collision term in the relativistic Boltzmann equa- 
tion to include nonlocal effects, and using Grad's 14-moment approx- 
imation for the single-particle distribution function, we derive evolu- 
tion equations for the relativistic dissipative fluid dynamics and com- 
£N| . pare them with the corresponding equations obtained in the standard 

Israel-Stewart and related approaches. Significance of this general- 
ization on hydrodynamic evolution is demonstrated in the framework 
of one-dimensional scaling expansion. 

O ! 1 Introduction 

Relativistic dissipative hydrodynamics has been quite successful in explain- 
ing the spectra and azimuthal anisotropy of particles produced in heavy-ion 
CNJ \ collisions at the RHIC and LHC [Tj. Apart from its applications, theoretical 

formulation of relativistic dissipative hydrodynamics is quite a challenging 
task by itself. The first-order dissipative fluid dynamics commonly known as 
relativistic Navier-Stokes (NS) theory involves parabolic differential equa- 
tions and suffers from acausality and instability. The second-order or Israel- 
Stewart (IS) theory [2] with its hyperbolic equations restores causality but 
\ may not guarantee stability. The correct formulation of relativistic dis- 

£f) ■ sipative fluid dynamics is still unresolved and is currently under intense 

investigation [3J3]. 

It is essential to note that all formulations of the second-order dissipa- 
tive hydrodynamics employing the Boltzmann equation (BE) make a strict 
assumption of locality in the configuration space for the collision term [2] . 
jj^ | This means, the collisions that increase or decrease the number of particles 

with a given momentum p, in an infinitesimal space-time volume element, 
are assumed to occur at the same point . This makes the collision inte- 
gral a purely local functional of the single-particle phase-space distribution 
function f(x,p) independent of the derivatives d^f. Although f(x,p) may 
not vary significantly over the length scale of a single collision event, its 
variation over the length scales extending over many inter-particle spacings 
may not be negligible. Including the gradients of f(x,p) in the collision 
term gives rise to different evolution equations for the dissipative quanti- 
ties. 

In this article, we provide a new formulation of the dissipative hydro- 
dynamic equations within kinetic theory by using a nonlocal collision term 
in the Boltzmann equation: p^d^f = C[f}. 
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INFINITESIMAL VOLUME ELEMENT IN A FLUID 




Figure 1: Collisions kk' — > pp' and pp' ~ > kk' separated by a spacetime 
interval £ M within an infinitesimal fluid element of size dR, containing a 
large number of particles represented by dots. 



2 Non-locality in collision term 

For two-body elastic collisions, the collision term is 

C[f] = X -l dp'dk dk' W pp ,^ kk ,(f k f k ,f p f p , - f p f p ,fhfk>), (1) 

where W pp >^ kk i is the collisional transition rate, f p = f(x,p), f p = l — 
rf{x,p) with r = 1, —1,0 for Fermi, Bose, and Boltzmann gas and dp = 
dp/[(2ir) 3 E p ]. The first and second terms in Eq. (JXJ) refer to the processes 
kk' — > pp' and pp' — > kk' , respectively. These processes are traditionally 
assumed to occur at the same space-time point x^ with an underlying 
assumption that f(x,p) is constant not only over a region characterizing a 
single collision, but also over an infinitesimal fluid element of size dR, large 
compared to the interparticle separation. It is important to note that the 
space-time points at which the above two processes occur may be separated 
by a small interval £ M within d 4 R (see Fig. [TJ. With this realistic viewpoint, 
the second term in Eq. (fjQ) involves f(x — £,,p)f(x — £,,p')f(x — £, k)f(x — 
£, k'), which on Taylor expansion up to second order in £ M , results in the 
modified BE [4] 

P^f = C[f] + B^A»f) + d^B^f). (2) 
The momentum dependence of coefficients, 

A" 

can be made explicit by expressing them in terms of the available tensors 
p^ and the metric g^ u = diag(l, -1, -1, -1) as A 1 = ap^ and B^ v = 
b\g^ v + bipVp 1 ' ' . The scalar coefficients a, b\ and 62 are functions of x^. 



Jdp'dkdk'^W pp ^ kk , fpl f k f k „ 

~ / dp'dk dk' e^W pp ,^ kk ,f p ,f k f k ,, (3) 



The conserved particle current, N^, and the energy- momentum tensor, 
T» v have the standard form g]: 



N» = j dp p»f = nu^ + ra", 

T" v = j dp p"p u f = eu^u v - (P + n)A^ + tt" u . (4) 



Conservation of current, d^N^ — and energy-momentum tensor, d t _ l T f * u — 
0, yield the fundamental evolution equations for n, e and it M 

Dn + nd^ + = 0, 

Dc + Ce + P + nja^-Tr^V^ti,,) = o, 
(e + P + n)L>u Q -V Q (P + n) + A^7r^ = 0. (5) 

Conservation of current and energy-momentum implies vanishing zeroth 
and first moments of the collision term C m \f\ leading to three constraint 
equations for the coefficients (0,61,62) 0], 

d^a = 0, d 2 (Moo) + dpdu (6 2 P"0 = 0, 

u a d^(b 2 I^ a ) + u a d 2 (b in u a ) = 0. (6) 

In order to obtain the second-order evolution equations for dissipative quan- 
tities, we consider second moment of the modified BE, Eq. © 

dp p a p fi p 1 d 1 ! =Jdp p a p? [C[f] + P W 7 (af) + d 2 (hf ) + (p ■ 9) 2 (6 2 / )] . 

(7) 

Using Grad's 14-moment approximation for the single particle distribution 
function in orthogonal basis, 

/ = /o + /o/o (A n n + X n n aP a + A ff7 r a/3 pV) , (8) 

we finally obtain the following evolution equations for the dissipative fluxes 

II = an NS - P^mtl + m n n ■ u - Ijmd ■ n - <5nnn6> + Ann" • Va 

+^n^^u^ u + A m it ■ u + Ai^ov^"" + (8 terms), (9) 
= an£ s - PnTnh^ + \ nn n v ^ - 5 nn n»6 + Z„nVTI - l^A^d^l 
-r„ n riu M - r„ 7r 7r MI/ u 1/ + \ nlx n v ^ iv + A nn nV + k n u^ v u v 
+A nu A^& y u 1v + (9 terms), (10) 
TT"" = OTT^ - faT*^ + T m n^u v > + Z™V<^> + K^lo^p 

+A^^ P + X l6 2 ^ + X2U ( ^b 2 + X3V^V>62, (11) 

in the usual notation 0]. The "8 terms" ("9 terms") involve second-order, 
linear scalar (vector) combinations of derivatives of 61, 6 2 . We observe that 
in the above equations, the nonlocal coefficients (a, 61, 6 2 ) modify the 
standard NS as well as the IS terms. This method is also able to generate 
all possible second-order terms allowed by symmetry. 
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Figure 2: Time evolution of temperature, shear pressure, inverse Reynolds 
number and parameters (61, 62) normalized to their initial values (left 
panel), and anisotropy parameter Pl/Pt (right panel). Initial values are 
T = 0.9 fm/c, T = 360 MeV, r)/s = 0.16, tt = 4t7/(3t ). 



3 Results and discussion 

To illustrate the numerical significance of the non-locality in the collision 
term through the new dissipative hydrodynamic equations derived here, 
we consider Bjorken evolution of a massless Boltzmann gas (e = 3P) at 
vanishing net baryon number density [5]. 

In terms of the coordinates (r, x,y, 77) where r = \Jt 2 — z 1 and 77 = 
tanh _1 (z/t), the initial four-velocity becomes = (1,0,0,0). In this sce- 
nario II = = tt/ 1 and the equation for tt = —T 2 -K m reduces to 

7r dir „ 4 , 7r , db2 , „, 

— + /?*— = ^— - A- -tM-A 12 
r w ar 3r t dr 

where the coefficients are 

A = a+ -7-5— , /3 n = -aP, V = ^-» A = 2a-2 — P. (13) 

ap/7 3 3ap?y 3ap7y 

The coupled differential equations ([5]), ((6]) and ()12|) are solved simultane- 
ously for a variety of initial conditions relevant for RHIC and LHC. 

Figure [5] (left panel) , shows the evolution of several quantities for a 
particular choice of initial conditions. T decreases monotonically to the 
crossover temperature T c ~ 170 MeV at time r ~ 10 fm/c. Parameters 
61, 62 vary smoothly to zero at large times indicating reduced but still 
significant presence of nonlocal effects at late times. This is also evident 
in Fig. [2] (right panel), where the pressure anisotropy Pl/Pt = (P — 
7r)/(P + 7r/2) shows marked deviation from IS. Although the shear pressure 
7r vanishes rapidly indicating approach to ideal fluid dynamics, the Pl/Pt 
is far from unity. Faster isotropization for initial a > may be attributed 
to a smaller effective shear viscosity in the modified NS equation. 
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Figure 3: Time evolution of Pl/Pt in IS (dashed lines) and the present 
work (solid lines) for isotropic initial pressure configuration. 

Figure |3] shows the evolution of Pl / Pr for isotropic initial pressure 
configuration, at various r\j s for the LHC energy regime. With small initial 
corrections (~ 10% to NS and ~ 20% to the IS terms) due to a, b\, 6 2 , 
nonlocal hydrodynamics (solid lines) exhibits appreciable deviation from 
the (local) IS theory (dashed lines). The above results clearly demonstrate 
the importance of the nonlocal effects, which should be incorporated in 
transport calculations as well. 

4 Summary 

To summarize, we have derived viscous hydrodynamic equations by intro- 
ducing a nonlocal generalization of the collision term in the Boltzmann 
equation. The Navier-Stokes as well as Israel-Stewart equations are modi- 
fied and new terms are obtained in the evolution equations of the dissipative 
quantities. The method presented is able to generate all possible terms that 
are allowed by symmetry. Within one-dimensional scaling expansion, we 
find that nonlocality of the collision term has a rather strong influence on 
the evolution of the viscous medium via hydrodynamic equations. 
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